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Abstract 


A literature review of recent research on two-phase flow in micro-channels is provided in this 
article. Researches on the micro-hydrodynamics concerned with two-phase gas—liquid adiabatic flow 
characteristics in both circular and non-circular micro-channels are discussed. This review aims to 
survey and identify new findings obtained from this attractive area, which may contribute to 
optimum design and process control of high performance miniature devices comprising extremely 
small channels. The results obtained from a number of previous studies show that the flow behaviors 
in the micro-channels deviate significantly from those in ordinarily sized channels. A clear 
understanding of the flow behaviors encountered is necessary to answer the fundamental questions 
related to the two-phase flow phenomena. Similar to the ordinarily sized channels, it is expected that 
the flow pattern, pressure drop and void fraction will influence the two-phase pressure drop, holdup, 
system stability, exchange rates of momentum, heat and mass during the phase-change heat transfer 
processes. It is obvious that publications available on micro-channels are relatively small compared 
to those for ordinarily sized channels and further systematic studies of two-phase gas—liquid 
adiabatic flow in micro-channels are required. 
© 2006 Elsevier Ltd. All rights reserved. 
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1. Introduction 


In the past years, studies on two-phase flow and heat transfer characteristics in 
micro-channel flow passages have become more necessary due to the rapid development of 
micro-scale devices used for several engineering applications including medical devices, high 
heat-flux compact heat exchangers, and cooling systems of various types of equipment such 
as high performance micro-electronics, supercomputers, high-powered lasers, and so on. 

In addition to the attractive phenomena caused by the significant role of the surface 
effects, this emerging field enables us to develop powerful miniature devices for operation 
in a domain which seemed to be unfeasible in the past. However, it is obvious that a 
comprehensive understanding is still lacking on the trends and parameters dominating the 
flow behavior in mini- and micro-channels. Thus, fundamental heat transfer and flow 
characteristics problems encountered in the development and processing of micro- 
electronic-mechanical systems (MEMS) have become a significant challenge for the design 
and control of these micro-systems. According to the requirement of new solutions for 
dissipating heat and keeping uniform temperature distributions in modern devices, a clear 
understanding of the major mechanisms affecting micro-scale heat transfer is essential for 
optimum design and process control of micro-systems. 

The two-phase flow and heat transfer characteristics in small channels such as micro- 
channels and mini channels are likely to be strongly dependent on the surface tension 
effects in addition to viscosity and inertia forces resulting in significant differences in the 
two-phase flow phenomena between ordinarily sized channels and small channels. Several 
investigators have proposed the criterion for a small channel based on different 
dimensionless parameters. 

Serizawa et al. [1] described the criterion for classification of micro-channel proposed by 
Suo and Griffith [2] as follows: 


A 
“> 
D >3.3, (1) 


where 4 is the Laplace constant defined by 


/ o 
ET 2 
i (PL — Po) ©) 


Brauner and Moalem-Maron [3] proposed the different criterion as follows: 


Eo < (2m0, (3) 
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Nomenclature 

C parameter dependent on Reynolds number 

Co distribution parameter 

D channel diameter, m 

Dg bubble diameter, m 

De equi-periphery diameter, m 

Da hydraulic diameter, m 

dP pressure drop, Pa 

dz unit length, m 

Eo Etovos number 

F time fraction 

‘fe friction factor 

G mass flux, kg/s m? 

g gravitational acceleration, m/s? 

J superficial velocity, m/s 

KEavg two-phase mixture’s average kinetic energy, kg/m s? 
length, m 

nı 1 for laminar flow; 0.25 for turbulent flow 

n 1 for laminar flow; 0.25 for turbulent flow 

Re Reynolds number 

s channel gap, m 

U average velocity, m/s 

V rising bubble velocity relative to the liquid phase, m/s 


We, | Weber number for vapor 

w channel width, m 

X parameter defined by Eq. (30) 
Xann parameter defined by Eq. (37) 
x flow quality 

a void fraction 

Am mean void fraction 

X Lockhart—Martinelli parameter 
h two-phase multiplier 

À Laplace constant 

u dynamic viscosity, kg/m s 

v kinematic viscosity, m?/s 

p density, kg/m? 

o surface tension, N/m 


Subscripts 

ann annular flow regime 

B bubble; slug bubble 

c contraction 

G gas phase, vapor phase 


GS gas slug 
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I gas-liquid interface 
int intermittent flow regime 
L liquid phase 


LO all-liquid 
LS liquid slug 


liq liquid flow regime 
s smaller channel 
T two-phase mixture 


UC unit cell 
VO all-vapor 
vap vapor flow regime 


where Eo is the Etovos number defined by 


— gL- po)D? 
=a. 


Eo (4) 

From Eqs. (1)-(4), o is surface tension, g is gravitational acceleration, D is channel 
diameter, py and pg are liquid and gas densities, respectively. 

Unfortunately, it is likely to be insufficient to identify the distinction between small 
channels and ordinarily sized channels by these equations. Some arbitrary channel 
classifications based on the hydraulic diameter D, have been proposed. Recently, 
Mehendale et al. [4] employed the hydraulic diameter as an important parameter for 
defining the heat exchangers as follows: 


Micro-heat exchanger: 
Meso-heat exchanger: 
Compact heat exchanger: 


1 um < Dp < 100 um 
100 um < D, < 1 mm 
1 mm < Dha <6 mm 


Conventional heat exchanger: D, > 6mm 


Based on the evolution of small flow channels used in engineering applications, however, 
Kandlikar [5] proposed the channel classification as follows: 


e Conventional channels: Dy>3mm 
e Mini-channels: 200 um < Dp < 3 mm 
e Micro-channels: 10 um < Dp < 200 um 


Although there have been a number of criteria based on the hydraulic diameter, the 
channel classifications available cannot relate the channel diameter to the fluid flow 
mechanisms. Further investigations should be conducted to meet a more general definition 
dealing with the channel classification. 

To obtain an appropriate design and process control of the various engineering systems 
employing these small channels, the two-phase flow characteristics as well as heat transfer 
mechanisms in various small channels, especially for micro-channels, should be clarified. In 
the past decade, there has been a relatively small amount of publications available for 
micro-channels compared to those for ordinarily sized channels. 
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The aim of this paper is to present a review of the most common types of research works 
done on two-phase gas-liquid adiabatic flow characteristics in micro-channels, and to 
propose the areas available for further investigations. 


2. Two-phase Gas-liquid adiabatic flow characteristics in micro-channels 


Up to now, there have been publications dealing with two-phase gas-liquid adiabatic 
flow characteristics in micro-channels with various geometrical configurations. However, 
the two most common types of such configurations employed by previous researchers are 
circular and non-circular micro-channels. 


2.1. Circular micro-channels 


Adiabatic two-phase air—water flow characteristics, including the two-phase flow pattern 
as well as the void fraction and two-phase frictional pressure drop, in circular micro- 
channels were experimentally studied by Triplett et al. [6,7]. For flow visualization results 
corresponding to their first paper, the flow patterns observed were bubbly, slug, churn, 
slug-annular and annular. These experimental results were compared with other 
experimental data obtained from previous investigators. The two-phase flow pattern 
transitions generally disagreed with the lines calculated from the existing models and 
correlations. 

According to their second paper, the measured values of void fraction and two-phase 
pressure drop in bubbly and slug flow regimes were found to be best predicted by the 
models based on the homogeneous mixture assumption. In the annular flow regime, 
however, the experimental void fractions and pressure drops were compared with the 
models as well as the available correlations, with poor agreement. 

Kawahara et al. [8] used a circular fused silica micro-channel with a diameter of 100 um 
to investigate two-phase flow characteristics. They observed many different flow patterns 
simultaneously in the circular channel under any given flow condition tested. Due to such 
simultaneous appearances, they used the probability of appearance of each flow pattern 
along with the time-averaged void fraction obtained from each flow condition to define the 
flow regimes as shown below. 


è Slug-ring flow was defined when the probability of “thin liquid film” was greater than 
that of “ring liquid film” and the time-averaged void fraction was less than 0.8. 

è Ring-slug flow corresponded to when the probability of “ring liquid film” was larger 
than that of “thin liquid film” and the time-averaged void fraction was less than 0.8. 

è Semi-annular flow was associated with flow in which the flow alternated mostly between 
“liquid alone” and “ring liquid film”. The time-averaged void fraction was greater than 
0.8. 

è Multiple flow included “liquid alone”, “thin liquid film”, “ring liquid film”, “thick 
liquid film” and “deformed interface”, and the time-averaged void fraction was less 
than 0.8. 

The flow patterns defined were used to develop a flow pattern map which was 
subsequently compared to other flow regime maps existing for two-phase flow of air—water 
in ~l mm diameter channels. It was noted that bubbly and churn flow patterns were 
absent in the experimental results due to the laminar nature of liquid flowing through the 
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micro-channel. The low value of the time-averaged void fraction, even at high gas flow 
rates, indicated the large slip ratios and weak momentum coupling between the phases. 
The friction factor in single-phase liquid flow and two-phase friction multiplier were also 
analyzed. The single-phase friction factor was found in reasonable agreement with the 
conventional correlation. The two-phase friction multiplier data correlated well with the 
Lockhart—Martinelli correlation based on a separated flow assumption, but the agreement 
between the data and homogeneous flow model was generally poor. 

Chen et al. [9] experimentally studied the characteristics of adiabatic two-phase flow of 
nitrogen—water in glass tubes with two different diameters, 1.0 and 1.5mm. Four common 
flow patterns including bubbly flow, slug flow, churn flow and annular flow were observed 
and a so-called bubble-train slug flow previously reported in a large channel under micro- 
gravity effects [10] was also observed and discussed in this work. The random manner of 
the number of bubbles counted in the bubble-train slug flow pattern was attributed to the 
turbulent nature. Moreover, increasing liquid superficial velocity tended to lower the 
number of bubbles. A modified drift flux model was found to be an appropriate approach 
to correlate bubble velocity for the three different flow patterns covering slug flow, bubble- 
train slug flow and churn flow. The correlation is expressed as 


Up = 0.932111, (5) 


where Upg is bubble velocity and jr is superficial velocity of the two-phase mixture. 
Furthermore, the experimental void fractions were well predicted by the use of this 
modified drift flux model in conjunction with the following fundamental equation: 


IG 
waa Q 
where « is void fraction and jg is gas superficial velocity. 

Serizawa et al. [1] investigated the visualization of the two-phase flow pattern in circular 
micro-channels. The flowing mixture of air and water in channels of 20, 25 and 100 um in 
diameter and that of steam and water in a channel of 50 um in diameter were conducted 
experimentally. Two-phase flow patterns obtained from both air—water and steam—water 
flows were quite similar and their detailed structures were described. The study confirmed 
that the surface wettability had a significant effect on the two-phase flow patterns in very 
small channels. It was also found that the air-water two-phase flow pattern transitions 
generally agreed with the lines predicted by the Mandhane’s correlation and the cross- 
sectional average void fraction concerned with bubbly and slug flows followed the values 
obtained from the Armand correlation [11] for an air—water system. 

Chung and Kawaji [12] performed an experiment in order to distinguish two-phase flow 
characteristics in micro-channels from those in mini-channels. Four different circular 
diameters ranging from 50 to 526um were employed to examine a scaling effect on 
nitrogen—water two-phase flow. The results including the flow patterns, void fraction and 
two-phase pressure drop were analyzed. 

For larger tubes tested (250 and 526 um), they observed flow patterns consistent with 
those obtained from mini-channels having diameters of around 1 mm. The corresponding 
flow patterns were bubbly, slug, churn, slug-annular and annular flow. The time-averaged 
void fraction agreed reasonably with the Armand-type correlation available for mini- 
channels. In the case of two-phase pressure drop, it was found that at a given 
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Lockhart—Martinelli parameter the two-phase friction multiplier tended to change with the 
variation of mass flux, similar to that appearing in mini-channels. 

For relatively small channels (49.5 and 100 um), only slug flow was observed under the 
experimental conditions tested. The void fraction and pressure drop characteristics were 
generally opposite to those of larger channels previously discussed. The relationship between 
the void fraction and volumetric quality was found to be non-linear and hence, the Armand- 
type correlation became unpredictable for these channel sizes. In addition, the frictional 
multiplier was independent of mass flux. Apparently, diameters between 100 and 250 um 
seemed to be in the range causing the notable change in two-phase flow mechanisms. 
Additionally, they modified Garimella et al.’s [13] model to predict the two-phase pressure 
drop for micro-channels having diameters of less than 100 um. In the proposed model, they 
assumed that the flow consisting of unit cells (UC) was completely intermittent as shown in 
Fig. 1. For each unit cell, there were two regions, including liquid single-phase flow region (or 
liquid slug), which had no entrained gas bubbles and the region occupied by the two-phase 
flow of liquid and gas. The latter region was formed by axially symmetrical gas slug 
surrounded by an annular liquid film with uniform thickness around the channel 
circumference. Because of the viscous effect, the velocity of the liquid film around the gas 
slug was assumed to be significantly low compared to that of liquid and gas slugs. 

In the single-phase region, the pressure gradient in the liquid slug is evaluated from the 
Darcy—Weisbach equation: 


dP _ p.Uis 
(Z) ahs 7D” (7) 


Fig. 1. Schematic diagram of unit cell for model of Chung and Kawaji [12]. 
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where (dP/dz);5 is the pressure gradient in the liquid slug, and fis is the friction factor for 
the liquid slug. The friction factor for fluid flow in a circular channel is given by 


64 
tis = Res > 
for laminar flow (Reis < 2100), (8) 
0.3164 


Tigo a5 
Rebe 
for turbulent flow (Reis < 100, 000) 


in a smooth channel. (9) 


Eq. (9) is known as Blasius equation. In the model, it was noted that, for the transition 
from laminar to turbulent flow (2100< Rers< 4000), the Blasius formula was also 
employed to obtain the friction factor. The Reynolds number of the liquid slug, Rer s, is 
defined as 


D 
Reys = PUP (10) 


HL 


where u; is the dynamic viscosity for liquid phase. In Eqs. (7) and (10), Ugs represents the 
average velocity of the liquid slug and is expressed by 


Uis =, (11) 
=% 


where ji is the superficial liquid velocity. 
In the two-phase region, the average velocity of the gas slug, Ugg, is taken to be equal to 
that of the gas phase, Ug 


Ugs = Ug (12) 
and 
Yeats (13) 
04 


The pressure gradient in the two-phase region, (dP/dz}r, is given by 


dP iat pg(Uas = Up? 
(=). PUGS pe 8 (14) 


where fos represents the friction factor for the gas slug and the determination of fgg is 
similar to that of fis. The Reynolds number of the gas slug, Regs, is given by 


pG(Ucs — U1)Dg 
HG í 


(15) 


Regs = 


where ug is the dynamic viscosity for gas phase. U; and Dg in Eqs. (14) and (15) represent 
the average velocity of the gas-liquid interface and bubble diameter, respectively. U; is 
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defined by Chung and Kawaji [12] as follows: 


(dP/dz) 


— T(p2 _ p2 
Ui = 5, P -DD (16) 


The bubble diameter is assumed to be 
Dz = 0.90D. (17) 


The unknown variables including Regs, Ur and (dP/dz); are obtained through 
Eqs. (14)-(16). Finally, the total frictional pressure gradient, (dP/dz)uc, across the whole 
unit cell is expressed by 


dP dP Les (=) Lis 
— = + , 18 
(i) UC (E) 7 Luc dz/15 Luc ne) 


where the relative lengths of the gas slug (Les/Luc) and liquid slug (Ly s/Luc) are 
approximated by 


Los ( D ) 

>= =u], 19 
Luc Dg (a) 
Lis Las 

— = 1 - —, 20 
Lyc Luc (20) 


Based on the micro-channels having diameters of 100 and 50 um, the experimental two- 
phase pressure drop data were predicted well by this slug flow model. 

Coleman and Krause [14] carried out an experiment to measure R-134a two-phase 
pressure drops due to contraction in circular micro-channel headers. The measured two- 
phase pressure drops were considerably high compared to some existing correlations. 
However, the authors suggested that the use of model developed by Schmidt and Friedel 
[15] in conjunction with the momentum density model was in reasonable agreement with 
the measured data. 

Two-phase air—water pressure changes in the circular channel, having sudden expansion 
and contraction sections, were experimentally studied by Abdelall et al. [16]. Two stainless 
steel tubes were connected to form such sections. The diameters of the larger tube and 
smaller tube were 1.6 and 0.84mm, respectively. The two-phase pressure changes were 
found to be considerably small when compared to those predicted by the homogeneous 
flow model. This indicated significant velocity slip in the neighborhood of the flow 
disturbance resulting from the area change. On the other hand, the pressure change 
data were well predicted by a simple one-dimensional (1D) flow theory with the slip ratio 
model based on the ideal annular flow regime with minimum entropy generation. They 
also developed an empirical correlation of two-phase multiplier for the flow area 
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contraction, dro <¢: 
Poe —0.7 
a = 120(y Retos) ~, (21) 


where Rezo, is the Reynolds number with the total flow in the liquid phase in smaller 
channel and the Lockhart—Martinelli parameter, y, is expressed by 


WA aA fa\™ 
x= (i) CS) GE)» o 
HG x PL 


where x is quality, wy and ug represent liquid and gas viscosities, respectively. Liquid and 
gas densities are denoted by py and pg, respectively. In Eq. (21), the subscript LO 
corresponds to the value with the total flow in the liquid phase and the subscript c and s 
correspond to contraction and the smaller channel, respectively. 


2.2. Non-circular micro-channels 


Triplett et al. [6,7] also conducted the same systematic experiments as discussed in the 
previous section to study gas-liquid two-phase flow characteristics in semi-triangular micro- 
channels. The results obtained were similar to those obtained from circular channels. 

Xu et al. [17] investigated an adiabatic co-current vertical two-phase flow of air and water in 
rectangular channels (12mm x 260mm) with mini-/micro-gaps of 0.3, 0.6 and 1.0mm. For 
channel gaps of 1.0 and 0.6mm, the flow patterns observed were similar to those found in 
ordinarily sized channels but peculiar flow structures which were different from typical two- 
phase flow patterns were observed in the channel having a micro-gap of 0.3mm. In the 0.3mm 
gap, they found that bubbly flow was never observed and new flow regimes including slug- 
droplet flow and annular-droplet flow were addressed. For the new flow regimes, it was observed 
that liquid droplets were always adhered to the wall surface and such droplets were blown by the 
moving gas phase. They also presented the models to predict the flow regime transitions in 
vertical rectangular channels. The models for bubbly flow to slug flow transition and slug flow to 
churn flow transition were developed based on Mishima and Ishii [18]. 

The transition from bubbly flow to slug flow is described by 


.  0.35./gDe 
1); (23) 


. 1 
Sr (az 


where C, is the distribution parameter and D, represents the equi-periphery diameter and 
is given by 
_ A4s+w) 


D. = : (24) 
T 


where s and w are the channel gap and width, respectively. In Eq. (23), the distribution 
parameter, Co, for rectangular channels is given by 


Co = 1.35 — 0.35, /PS. (25) 
PL 


The slug-churn flow transition took place when 


> Om. (26) 
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The void fraction, «, is given by 


ja 
a=- (27) 
Cola +j) + V 


where V is the rising bubble velocity relative to the liquid phase and is expressed by 


V =0.35,/gDe. (28) 
The slug bubble mean void fraction &m, shown in Eq. (26), is defined as 
én = 1—0.813x"™, (29) 
where 


_ PL : ; 
X = Jaar (Ce - Mic +) + 0.35V/9D;]. (30) 


Lpg and Ap are the slug bubble length and density difference between the two phases, 
respectively. 

Based on a void fraction of 0.75, the transition of the annular flow took place when gas 
and liquid superficial velocities satisfied the following relationship: 


1.3374 — 4, 
(P — pa)gDn + 12861 (424) “pu, 


9.24C6 (2) "pg i +1.5(2) ”] Sa (31) 


where Cy and Cg are dependent on the Reynolds number. 

In Eq. (31), v_ and vg, are kinematic liquid and gas viscosities, respectively. nı or nz is 
equal to unity for a laminar flow and nı or m is equal to 0.25 for a turbulent flow. 
However, the models developed were found to be predictable for gaps larger than 0.6 mm. 
For the 0.3mm micro-gap, they suggested that an appropriate theory consistent with 
peculiar flow regimes observed was needed to develop the flow regime transition models. 

Zhao and Bi [19,20] performed experiments to study two-phase flow patterns in vertical 
upward air—water flow in equilateral triangular channels having hydraulic diameters of 
2.886, 1.443 and 0.866mm. The images of the flow patterns were taken by a high-speed 
motion analyzer along with the transient pressure-drop measurements. The flow patterns, 
including dispersed bubbly flow, slug flow, churn flow and annular flow, were observed in 
the two larger triangular channels. These observations were similar to those generally 
encountered in the ordinarily sized tubes placed in vertical direction. However, it was 
interesting to note that a new flow pattern, namely capillary bubbly flow, was identified in 
the smallest channel in addition to the typical flow patterns, excluding the dispersed bubbly 
flow. Another interesting point to note was that the slug-churn flow and churn—annular 
flow transition lines shifted to the right as the hydraulic diameter decreased. In addition, 
the previous flow regime transition models being available for upward gas-liquid two- 
phase flow in the ordinarily sized tubes could not be used to predict the flow transitions. 

According to their second paper, the experimental investigation for analyzing gas 
velocity, void fraction, and pressure drop was carried out in the same experimental setup. 
They reported that, based on the distribution of the parameter suggested by Ishii [21] and 
based on negligible drift velocity, the drift-flux model was in reasonable agreement with the 
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experimental data for gas velocity obtained from all three test sections. In addition, they 
obtained a correlation for predicting the void fraction and single-phase pressure drop. The 
void fraction correlation developed was similar to that suggested by Armand and Treschev 
[11]. It was found that the Lockhart—Martinelli correlation was in good agreement with the 
two-phase pressure drop data when the proposed new single-phase friction factor 
correlation was employed. 

Chung et al. [22] explored nitrogen gas—water two-phase flow characteristics in a square 
channel having a hydraulic diameter of 96 um and compared the results with those for a 
100 um circular channel employed by Kawahara et al. [8]. Similar results including pressure 
drop and void fraction data were obtained, regardless of the channel geometry. In both 
channels, the Lockhart—Martinelli correlation appeared to do rather well when compared with 
the pressure drop data while the void fraction data were found to be in poor agreement with 
the homogeneous flow model as well as the Armand-type correlation. They also discussed the 
effect of channel shape on the transition boundaries for flow patterns observed. 

A flow visualization study to clarify the flow patterns of vertical upward gas-liquid two- 
phase flow in rectangular mini-channels with hydraulic diameter ranging from 1.95 to 
5.58mm was carried out by Satitchaicharoen and Wongwises [23]. Air—water, air—20 wt% 
glycerol solution, and air-40 wt% glycerol solution were used as working fluids. In the 
experiments, they employed various rectangular test sections: 20mm x 2mm, 
40mm x 1 mm, 40mm x 2mm, 40mm x 3mm and 60mm x 2mm at an equal length of 
1m. The flow phenomena, which were bubbly flow, cap-bubbly flow, slug flow, churn flow 
and annular flow, were observed and recorded by a high-speed camera. The effects of gap 
size, channel width and liquid viscosity on the flow pattern transitions were discussed. 

English and Kandlikar [24] carried out experiments to explore the effect of surfactants 
on two-phase pressure drop of air—water adiabatic flow in a 1 mm square channel. The 
conditions they tested were based on the operating control in a proton exchange membrane 
(PEM) fuel cell. The surfactant Triton DF-12 was employed to reduce the surface 
tension to 0.034N/m. The application of the surfactant solutions appeared to have no 
effect on the pressure drop due to many possible reasons, such as characteristics of the flow 
pattern affecting interaction between the two phases, insufficient surface tension for an 
observable change, and temperature variation. The accuracy of existing models for two- 
phase pressure drop prediction was evaluated. Based on the Mishima and Hibiki model 
[25], English and Kandlikar [24] also developed a new model applicable to low mass fluxes, 
high mass qualities and annular flow by applying the Chisholm value for laminar—laminar 
flow of 5 instead of the original constant of 21. 

Two-phase flow characteristics including flow pattern and pressure drop of ethanol-CO, 
flowing through a converging or diverging rectangular micro-channel were explored by 
Hwang et al. [26]. Different from previous works associated with the flow in micro- 
channels, they observed stratified wavy flows at the entrance of the converging micro- 
channel. Peculiar flow patterns and bubble interactions, which were due to the presence of 
pressure and acceleration effects in the variable cross-section channels, were also discussed. 
For the pressure drop, they found that the channel pressure drop in the diverging micro- 
channel was substantially small in comparison to that in the converging micro-channel. 
However, the two-phase frictional multiplier plotted against the Lockhart—Martinelli 
parameter for both channels was independent of the liquid flow rate. This result was in 
agreement with that observed in nitrogen—water two-phase flow in circular micro-channels 
investigated by Chung and Kawaji [12]. 
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Due to different flow regimes raised in the same micro-channel under a given flow 
condition, Jassim and Newell [27] proposed void fraction and pressure drop models based 
on the probabilistic flow regime map. This modeling technique was developed in such a 
way that flow visualization data were incorporated into the models developed for 
predicting pressure drop as well as the void fraction for air—water at 20 °C and refrigerants 
including R410A and R134a at 10°C with mass flux values ranging between 50 and 300 kg/ 
m?s. A test section with 6-port micro-channels of 1.54mm in hydraulic diameter was used 
to obtain time fraction data for different flow regimes including liquid, intermittent, 
annular, and vapor flows. The proposed probabilistic model for predicting total pressure 
drop in micro-channels can be expressed with many flow regime time fraction functions 
and is given by 


di gp APY. ge (EP) ge HEP) ae PE oa 
dz total =a dz liq mi dz int es dz vap a dz ii 


where Fig, Fins Fvap and Fann represent time fraction curve fits in liquid, intermittent, 
vapor and annular flow regimes, respectively. The pressure drop for single-phase flow, 
including that for liquid flow (dP/dz)iq and that for vapor flow (dP/dz)yap, is determined 
from the conventional relation: 


dP Gl 
dP = aN 33 
ah or vap Ta (5) = 


where G is mass flux. The pressure drop model for intermittent flow assuming a 
homogeneous density, pr, and employing the average kinetic energy, KE,yo, is given by 


dP 1 
— |) =0.045—KEaye, 
(h = aD e 


h 
where 
2 
E 
2pr 
and 
x lee ae 
PT = = + ( ? ; (34) 
PG PL 
The annular flow pressure drop model is evaluated as follows: 
dP > (dP 
a ann 7 Pre (i) vo. Pa 


where Po is two-phase multiplier and is expressed by 


Pro = exp (—0.046X ann) + 0.22[exp (—0.002X ann) — exp (—7X ann)], (36) 


| ra 0.9 m 
(+ wa) (Fe) | e» 
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where We, in Eq. (37) represents the Weber number for vapor and is defined as 


2 
We, = LEST Pol (38) 


In Eq. (35), the vapor only pressure drop is given by 


dP 1 æ 
(£) = 2fvo eee (39) 
Z/ vo h PG 
The probabilistic void fraction model proposed is also expressed with many flow regime 
time fraction functions as shown by 
“total = FrigK%iq + FintOint + Fyap%vap + Fann%anns (40) 


where oq represents the void fraction for liquid only flow and is equal to 0, while &vap 
represents the vapor only flow void fraction which is equal to 1. The intermittent flow void 
fraction, čin and annular flow void fraction, «ann, are defined as follows: 


(0.833 + 0.167x)x(1/pg) 
(1 -= x)(1/pL) + x(1/ec) ’ 


j 0.97 —0-06 
; PL 
HG n 


Pressure drop and void fraction predictions based on the probabilistic flow regime map 
were found to be in accordance with the measured data obtained in micro-channels. 
However, care should be taken when the models are needed to be applied under conditions 
and geometries other than those in which the probabilistic flow regime maps were 
established. 


(41) 


Xint = 


Cann = 


3. Conclusions 


As discussed above, this review reveals that there is only a small amount of publications 
available for the two-phase gas-liquid adiabatic flow in micro-channels. The existing 
literature is still lacking in various aspects, including geometry of mixing chamber 
configuration, fluid properties, channel orientation and relevant flow models. 

It is obvious that the status of current researches on phase-change heat transfer and two- 
phase flow phenomena in micro-channels is still in infant stage and a great deal of 
systematic investigations will become substantially more necessary to meet general 
conclusions needed for the appropriate design and process control of several engineering 
applications. 
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